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We compute the gluon and ghost propagators of Yang-Mills theory in linear covariant gauges from 
the coupled system of Dyson-Schwinger equations. For small values of the gauge fixing parameter 
(, ~ 0.1, the deviations to the Landau gauge already become clearly visible. For the ghost dressing 
function, this is reflected in a logarithmic infrared suppression. Also, the gluon dressing function 
changes - but only quantitatively - and the gluon propagator remains finite at zero momentum. 
From the results, a running coupling is extracted. 
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I. INTRODUCTION 

Quantum chromodynamics is a theory rich with physical 
phenomena, e.g., asymptotic freedom, confinement, dy¬ 
namical mass creation, to name only a few well-known ex¬ 
amples. The investigation of its low energy regime faces 
the challenge of being nonperturbative and requires cor¬ 
responding methods like lattice simulations or functional 
methods. 

In this paper we will investigate Yang-Mills theory - 
i.e., we neglect all quark effects - with Dyson-Schwinger 
equations (DSEs). The basic quantities in this approach 
are correlation functions of helds. The lowest ones, the 
propagators, are particularly important, as they encode 
most of the nonperturbative dynamics. As correlation 
functions are gauge dependent quantities, the choice of 
a gauge is required. This allows, to some extent, to 
avoid or alleviate some difficulties. For example, differ¬ 
ent gauges feature a different number of fields and/or 
interactions. Thus gauges with a minimal number of 
both have a technical advantage over others. Indeed, 
the gauge which was most prominently used during the 
last decades fulhlls this criterion: The Landau gauge has 
two fields (gluon and ghost) and three primitively diver¬ 
gent vertices. Calculations in this gauge have advanced 
as far as the self-consistent calculation of its primitively 
divergent vertices m- Of course there are also other 
reasons for the attractiveness of the Landau gauge. One 
is the fact that the transverse correlation functions form 
a closed system and decouple from the longitudinal ones 
[5] . Historically the simple form of the ghost-gluon vertex 
played an important role as it provided the entry point 
to the coupled system of gluon and ghost propagators 
[^. Since then, much work on the propagators has been 
done; see, for example, [HEIEHIT]. 

Other gauges besides the Landau gauge have also been 
used in the past. Among them, the Coulomb gauge is 
the gauge where the study of the elementary Green func¬ 
tions has also progressed rather far [IHH2I]- Another 
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example is the maximal Abelian gauge, the correlation 
functions of which have also been studied to some ex¬ 
tent, e.g., [m [55]. And, last but not least, there are 
the linear covariant gauges, the end point of which is the 
Landau gauge. Within the functional framework they 
have been investigated only occasionally, e.g., [231 155] . 
Results from different gauges offer the possibility of test¬ 
ing the gauge (in)dependence of observables. Naturally, 
observables are gauge independent, but truncating the 
underlying functional equations can spoil this property. 
In this respect linear covariant gauges play a special role, 
as they allow us to change the gauge continuously and 
have the Landau gauge as their end point, which is well 
studied. 

At first sight, the extension to a nonvanishing value of the 
gauge fixing parameter of the linear covariant gauges, 
may seem rather straightforward once the necessary tech¬ 
niques are mastered. However, the little knowledge we 
have about this family of gauges for providing a guide in 
the construction of a truncation for functional equations 
is mostly based on an extrapolation from the Landau 
gauge. In addition, the longitudinal parts of the correla¬ 
tion functions have to be considered. 

In this context it is interesting to note that with lattice 
calculations the extension to nonzero ^ is also a nontriv¬ 
ial issue. However, the main reason for this is the fact 
that conventional gauge fixing techniques are based on 
the extremization of a functional; see Eq. Q below. 
Such functionals exist for the Landau gauge, the maxi¬ 
mal Abelian gauge, and the Coulomb gauge, but for lin¬ 
ear covariant gauges one can show that such a functional 
cannot be constructed in the conventional way and that 
alternative methods are required |55H55] ; see Ref. [50] for 
an overview on this topic. 

Here we present an exploratory investigation of the prop¬ 
agators in linear covariant gauges using Dyson-Schwinger 
equations. The truncation we use is motivated by its 
successful application in the Landau gauge. Lacking 
information on the transverse parts of the vertices we 
construct models that smoothly connect to the Landau 
gauge but contain only the minimally required ingredi¬ 
ents like the correct ultraviolet (UV) behavior. For the 
longitudinal parts we rely on information provided by 
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Ward identities. The resulting system of equations is 
then solved self-consistently. 

The article is organized as follows. Yang-Mills theory in 
linear covariant gauges is discussed in Sec.|nj The Dyson- 
Schwinger equations for the propagators are introduced 
in Sec. |III[ where the IR behavior of the ghost propagator 
is also investigated. The numerical results are presented 
in Sec. |IV| and in Sec.|V]we summarize. In the Appendix 
some details of the Dyson-Schwinger equations are dis¬ 
cussed. 


II. LINEAR COVARIANT GAUGES 

The action of Yang-Mills theory, S'ym = j / 

is invariant under gauge transformations of the gauge 

field 

A“ ^ Df = S‘^% + g (1) 

The Landau gauge is defined as minimizing the norm of 
the gauge fields 

= ( 2 ) 

with respect to gauge transformations. Any gauge 
field configuration obeying 9^^^“ = 0 fulfills this crite¬ 
rion. However, nonperturbatively this gauge fixing is not 
unique and the issue of Gribov copies arises [SI]. The 
restriction to the hyperplane = 0 introduces aux¬ 

iliary fields called ghosts via localization of the Jacobian 
determinant: 

= I d^x . (3) 

In the limit ^ = 0 the condition = 0 is strictly 

enforced. This is the case of the Landau gauge. ^ > 0 
corresponds to a Gaussian distribution along the gauge 
orbit centered at = 0 with a width of Gri¬ 

bov copies manifest then in the existence of several such 
Gaussians. 

The dressed gluon propagator is parametrized as 

DfAp)=Df/{p) + ^5‘^^^-^ (4) 

with 

fiab 

= -^Z{p^)P,Ap) ( 5 ) 

being the transverse part. P^v{p) = g^ll/ — P^iPv/p^ is 
the transverse projector. As is known from the Slavnov- 
Taylor identity for the gluon propagator, the longitudi¬ 
nal part stays bare. In perturbative calculations, where 
Z{p^) = I is used, the Feynman gauge with ^ = 1 sim¬ 
plifies the calculations, but in a nonperturbative setting 
the dressing of the transverse part negates that and the 


Green function 

Anomalous dimension 

Ghost propagator 

A - _9z:M 

^ 44 

Gluon propagator 

1 

1 

II 

c- 

Ghost-gluon vertex 

_ 3^ 

Three-gluon vertex 

...Sg _ 17-9« 

1 44 


TABLE I. The anomalous dimensions of the propagators and 
vertices. 


Feynman gauge is not simpler. Another peculiar value 
for ^ is 3, for which the ghost self-energy becomes finite; 
see Tab. |Tj This is the Yennie gauge. The unitary gauge 
is ^ oo. 

The inverse of the propagator, the two-point function, is 
written as 

r“t(T) = (rf/(p) + r'<5“ W) (6) 

with 

rf’^ip)=6‘^^Z-^{p^)P,Up)p^. (7) 

The ghost propagator D‘^’’{p) is parametrized by 

= ( 8 ) 

The two dressing functions of the ghost-gluon vertex are 
taken as the parts transverse and longitudinal with re¬ 
spect to the gluon leg: 

{Df{k‘^]p'^,q^)Pf,^{k)pu + Df{k‘^-,p^. 

(9) 

The momenta k, p, q are those of the gluon, the 
antighost and the ghost. The bare vertex T‘^^'^’^^\k-,p,q) 
is obtained by setting D^{k"^]p^ ,q^) = 1 and 

,q^) =p-k/k‘^. In the Landau gauge the ghost- 
gluon vertex is often used without radiative corrections. 
In particular, it is known that it is UV finite |32j . This is 
no longer true for ^ > 0 and we implement the running 
of the vertex via the following model: 

Df{k^-p\q^) = G{p^ri^Z(fYi\ ( 10 ) 

p^ is -\- q^ Pk"^)/2. The exponents of® and /3f® can be 
found in Tab. [T^ in the Appendix where their calculation 
is also discussed since they are the same ones used for 
the renormalization group improvement. Using Tab.[l] it 
can be easily checked that this expression has the correct 
UV behavior. 

The full three-gluon vertex has 14 tensors, of which four 
form the transverse subspace. Motivated by results in 
the Landau gauge, where it was shown by explicit calcu¬ 
lations that all dressing functions of the transverse ten¬ 
sors are negligible except for the contribution from the 
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FIG. 1. Truncated two-point Dyson-Schwinger equations. All 
internal propagators are dressed. Thick blobs denote dressed 
vertices. Wiggly lines are gluons, dashed ones ghosts. Plots 
were created with Jaxodraw m- 


tree-level tensor [5], we restrict ourselves to this tensor 
given by 

tgr^mr-q)^g,, + peTm.). ( 11 ) 

The transverse part of the full three-gluon vertex is then 
written as the transversely projected bare tensor with a 
dressing ,p^,q^) given by 




For the longitudinal parts of the vertices we use Ward 
identities. They are obtained from the invariance of the 
path integral under gauge transformations, see Refs. [31- 


1311 fo'' details: 



is the Ward operator. From this identity we calculate the 
Ward identities for the three-point functions by applying 
further derivatives. As an approximation we keep only 
the terms of lowest order in g. The resulting expressions 
are then used for the longitudinally projected vertices in 
the DSEs. For the ghost-gluon vertex we obtain 

(15) 

and for the three-gluon vertex 

=5/“''^ - ^lp{q)) ■ (16) 

III. THE PROPAGATOR DYSON-SCHWINGER 
EQUATIONS 

Based on the invariance of the path integral under trans¬ 
lations of the fields, the equations of motion for all cor¬ 
relation functions can be derived; see, e.g., 0 [38] for 
details. On a formal level the DSEs in linear covariant 
gauges look the same as in the Landau gauge, viz.. 


(D“'’(p)) ^ = -Z,p^+g^Z, D^^{p + q)D;i{q)rl‘‘\-q;p + q-p)Tl^^’^°\q;p,-q-p) (17) 

Jq 

+ r^S‘^%p, + / DYq)D^Yp + qK^^^i-P^pP q, -g)r““’(°)(T; q,- p- q) 


- y ^ DfAq)D%{p + q)TYi{-p,p + 


see also Fig. stands for / d"‘g/(27r)'^ and Zi and 

Zi are the renormalization constants of the ghost-gluon 
and three-gluon vertices, respectively, whereas Z^ is the 
renormalization constant of the ghost propagator and Z^ 
that for the gluon propagator. The dots represent two- 
loop terms and the tadpole diagram. The latter is typ¬ 
ically neglected as it does not contribute perturbatively 
and its nonperturbative impact is small at best [39) . The 


) -q)Y7a°^ (t; 9, -T - 9) + • ■ •, (18) 


two-loop terms, on the other hand, are important in 
the midmomentum regime, at least in the Landau gauge 
0130 ]. However, for the moment we leave these contribu¬ 
tions aside. Before we can aim at quantitative results we 
would need more detailed information about the vertices 
in general, which is not available. 

We bring the DSEs now into a form amenable to numer¬ 
ical computation. We start with the ghost propagator 
DSE ( [T7| ). Plugging in explicit expressions for the prop¬ 
agators and vertices, it becomes 
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G{x)-^ =Z3 + 


g^N.Zi f G{z) 


IqXy'^z 


{{x^ + {y- zf - 2x{y + z))Z{y)Df{y, z, x) - 2^y{x + y- z)Df{y, z, x)) . (19) 


Here and in the following, we use x = p^, y = and z = {p + q)^. Switching to hyperspherical coordinates, we can 
integrate out two angles and the integrand reduces to two terms: 


N-i _ 9 [, , gsin((/3)2 / s\n{ipfD^^<^{y,z,x) ^qcosp 


G{x) — Z^ + 


47^3 


dq d(p- 


ZA{y) 




p 


Dr^iy,z,x)]. 


zZciz) 

The final expression is obtained by plugging in the Ward identity for the ghost-gluon vertex from Eq. (15): 

g^N.Z, 


( 20 ) 


G (x) — Z^ + 


87r3 


dy ^^ 3in(v;) ^ jjAcc ^ x)Z{y) + i ^ cos(pf ), (21) 

z \ qG{x) 


where one term vanished due to the angle integration. 


From this expression the infrared (IR) behavior of the 
ghost propagator can be inferred. As we will argue, the 
term proportional to ^ is IR leading, which entails a dif¬ 
ferent IR behavior from the Landau gauge. In the Lan¬ 
dau gauge, only the first term appears, which goes to 
zero for p —)• 0 since the product G{z)Z{y)D^^‘^{y,z,x) 
vanishes faster than the rest of the integrand. If that 
were different for ^ > 0, the qualitative IR behavior of 
the propagators or the ghost-gluon vertex would need to 
change drastically. We will assume that this is not the 
case and will support this by an explicit calculation. For 
low external momentum p, the third part becomes 




g ^N,Z, 

8tt^ 


dydip ^"^^^^^ cos(p)^G(2/). 

y 


( 22 ) 


If G{y) were IR divergent or constant this would lead 
to an IR divergence in the integral. In turn this would 
mean that G{x) on the left-hand side would vanish in 
contradiction to the original assumption. Thus the ghost 
dressing function must be IR suppressed to make the 
integrand convergent. The numerical results will indeed 
show this behavior. 

We now turn to the gluon propagator DSE. To transform 
it into a scalar equation we project it with the transverse 
projector. Note that then only the transverse part of the 
ghost-gluon vertex appears. Splitting the gluon loop by 
orders in the DSE reads 

Z 3 + g^ N, Zi f G(j/)G(z)i^f (x; y, z)Df{x; y, z) 

Jq 

+ g^ N, ^ (z(y)Z(z)if|'(cr; p, z)D^^{x, y, z) 

+ Z{x)-^'^ . (23) 


The kernels are given in Eq. (26) in the Appendix. Note 


that in the gluon dressing function appears in var¬ 

ious combinations. 


I 

Before the system of equations can be solved some further 
modifications are required. They are related to the UV 
behavior of the equations, which should consistently pro¬ 
duce the one-loop resummed perturbation theory. How¬ 
ever, the truncation spoils this property. Since the em¬ 
ployed Newton method is sensitive to a consistent UV 
behavior, we restore this property by adding renormal¬ 
ization group improvement terms [T]. The procedure on 
how to do this is described in the Appendix. The issue 
of spurious UV divergences which are subtracted pertur- 
batively is also discussed there [39] . The final equations 
including these modifications can be found in the Ap¬ 
pendix in Eqs. (28) and (29). 


An interesting question with regard to the IR behavior 
of Green functions is if a scaling solution exists also for 
nonzero Such a solution is characterized by power laws 
for the dressing functions ESI]. A DSE analysis using a 
bare vertex approximation found that this is not possible 
if one requires that the longitudinal gluon propagator 
stays bare [Mj. This also holds if dressed vertices are 
considered [42] . The only way out would be nontrivial 
cancelations in the loop diagrams, which, however, we 
do not find here. However, given the analytic result that 
the ghost dressing function is IR suppressed, it is hard 
to see how a scaling solution could be realized. 


IV. RESULTS 

Using the framework provided by DoFun 1381 143] and 
CrasyDSE [33] the system of the DSEs (28) and (29) 


is solved using a standard Newton method, see, e.g., [3 
H51H5] . Solutions for various values of ^ including 0 are 
determined. It turned out that for a nonzero gauge fixing 
parameter, a higher precision is required and that with 
the currently employed setup the value of ^ is limited to 
^ < 0.4. For larger values, numerical artifacts become 
too large. 

Our results are shown in Fig. [^ Even though Yang- 
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FIG. 2. Ghost (left) and gluon (right) dressing functions for various values of 


a99(p2) 



a^V) 



p[GeV] 


FIG. 3. Running couplings from the ghost-gluon (left panel) and the three-gluon (right panel) vertices. 


Mills theory does not have a physical scale, one can use 
the string tension or quark potential methods to set one 
artificially. To inherit this scale, we use the distinct bump 
in the Landau gauge gluon dressing function as a marker. 
Using truncation schemes including the vertices one can 
obtain with this method very good agreement with the 
running coupling in the universal perturbative regime [T] . 
Here we used lattice results from Ref. [IB] and calculated 
the scale setting factor from the Landau gauge result. It 
was then applied for all values of 

First differences to the Landau gauge are seen in the 
ghost dressing function which starts to bend down in the 
IR. This effect, which can be seen already for very small 
values of the gauge fixing parameter, ^ ~ 0.001, where 
the gluon dressing function is still unaffected, is in agree¬ 
ment with the qualitative IR analysis in Sec. m The 
IR suppression is logarithmic and not as strong as 
which would make the ghost propagator IR finite. This 
behavior was found in an earlier DSE analysis |25j . In an¬ 
other study using a variational method the ghost dressing 
function was found to be constant at low momenta m- 
Around ^ = 0.1 the gluon dressing also starts to change 
and the bump in the midmomentum regime becomes big¬ 
ger. At this value of the gauge fixing parameter, the ghost 
dressing function deviates in the IR already severely from 
the Landau gauge. Raising ^ further, the ghost dressing 


function becomes even more IR suppressed. The effect in 
the gluon dressing is most drastic in the midmomentum 
regime, where the bump gets enlarged. At small mo¬ 
menta the gluon dressing function always vanishes like 
so that the propagator is IR finite. This is in agree¬ 
ment with the result of Ref. m- The UV behavior for 
both propagators changes as expected. 

While for small values of ^ the proximity to the well- 
known Landau gauge most likely makes the results trust¬ 
worthy, it is unclear how well the employed approxima¬ 
tions work for larger values of In particular, the shift 
of the bump in the gluon dressing function to higher mo¬ 
menta and the large increase in height indicates that the 
employed truncation may still lack some important fea¬ 
tures. An obvious possibility for improving the trunca¬ 
tion is to use more realistic vertex dressings. The em¬ 
ployed models are rather simple and do not capture any 
nontrivial structures in the nonperturbative regime. 

We also calculate the running couplings as extracted from 
the ghost-gluon and three-gluon vertices (B] |3H] : 

a8S(p2) = a(p2) G^{p^)Z{p^), (24) 

a3S(p2) = q,(^2) ^ 25 ) 

Note that the ghost-gluon vertex enters explicitly, in con¬ 
trast to the Landau gauge, where it is UV finite. The 
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results are shown in Fig. In the perturbative regime 
they all agree as expected from universality. The large 
bump in the gluon dressing function leads to a rise at 
intermediate momenta for a larger For increasing val¬ 
ues of the bump in the coupling also moves to higher 
values. Unfortunately, our results do not show any sign 
of a slowing down of that movement, so that for large 
enough values of ^ this will be in conflict with pertur¬ 
bative universality. This hints again at a shortcoming of 
the vertex models for larger values of 
Finally, we note that a desired property of our solution 
is that it fulfills a confinement criterion based on the 
Polyakov loop potential gaED]. The important property 
is that the gluon propagator is IR suppressed relative to 
the ghost propagator, which then leads to a confining 
Polyakov loop potential. 


Thus the ghost propagator has a different IR behavior 
from the Landau gauge. The gluon propagator, on the 
other hand, qualitatively has the same IR behavior; viz., 
it becomes constant. With an increasing value of the 
gauge fixing parameter the bump in the gluon dress¬ 
ing function shifts towards larger momenta and becomes 
higher. However, this may be a truncation artifact. To 
investigate this issue further dedicated investigations of 
the vertices - or at least the use of improved models - 
would be required. 

Note added: Recently Ref. m appeared, where the prop¬ 
agators of linear covariant gauges were investigated. The 
numerical results obtained for the ghost propagator are 
in qualitative agreement with ours. 
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V. SUMMARY AND CONCLUSIONS 

In this work we presented results for the propagators of 
Yang-Mills theory in linear covariant gauges from Dyson- 
Schwinger equations. The truncation was kept quite sim¬ 
ple, with the only dynamic quantities being the propaga¬ 
tors. For the vertices, models were used for the transverse 
parts that respect the correct UV behavior and the lon¬ 
gitudinal parts were taken from their Ward identities in 
leading order in g. From the results for the propagators 
a running coupling was extracted. 

Analytically, we could show from the ghost propagator 
DSE that the ghost dressing function must be IR sup¬ 
pressed. The numerical results confirmed that behavior. 
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Appendix: Details on the Dyson-Schwinger 
equations of the propagators 


The kernels for the gluon propagator DSE in Eq. 
are given by 


(23) 


K^z^{x-,y,z) 

Kt{x-,y,z) 

Kp^{x;y,z) 


- 2x{y + z) + {y- z^) 

12x‘^yz 

{x^ — 2x{y + z) + {y — z)^) (x^ -I- 10x(2/ + z) + y^ + lOyz + z^) 

24a;2y2z2 ’ 

(x^(y -I- z) -I- x^ (9?/^ - Ayz + 9z^) -|- x {-9y^ + y'^z + yz^ - 9z^) - {y - z)^ -I- z^)) 

24x2y2z2 

Z(z)(x - z) (x^ - 2x{y - 5z) + {y - zf) ^ Z{y){x - y) (x^ -b 2x(5y - z) + {y - zf) 
2Axy‘^ z"^ Z {x) 2Axy'^ z'^ Z (x) 

{x^ - 2x{y + z) + {y - zf) 

2Ay^z^^''^" 


(26a) 

(26b) 


(26c) 

(26d) 


To solve the system of equations, a standard New¬ 
ton method is employed. A stable iteration is only 
achieved if the anomalous dimensions are reproduced 
self-consistently. However, one-loop resummed pertur¬ 
bation theory requires a resummation of the diagrams 
not included in the truncation. In the Landau gauge, 
one possible remedy is to modify the vertex models such 


that the correct running is produced, or, in other words, 
to use momentum dependent renormalization constants 
[I1IS1I3S]. Because of the presence of the mixed terms 
in the kernel this procedure needs to be modified: 

Every term gets multiplied by a so-called renormalization 
group improvement factor to obtain the correct one-loop 
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running [T]. This factor is generically 

F{a,(3-,p^) = G{p^rZ{py, (27) 

where p^ is {x + y + z)/2. This choice ensures that for 
large loop momenta q the argument becomes y while be¬ 
ing symmetric in all momenta. The exponents a and 
P are determined such that the logarithmic running in 
the UV is correct. In addition we would require, as in 


the Landau gauge, that F(a,/3;p^) would become con¬ 
stant in the IR if the ghost propagator were IR constant. 
However, since the ghost propagator is not IR constant 
the latter condition should be improved in future calcula¬ 
tions, particularly as it can have a quantitative influence 
on the results [a [39]. However, we adopt it here due to 
its simplicity. The final DSEs then read 


G-^{x) = Z3 + 




dy dip 


sin((/?)^ 


- sin{p)‘^z, x)Z{y)F{af, 


+ ^F(af,/3f;p-2) (^^^+cos((^)2j jG(,), 


Z-\p^) = Z^+g^N, [ G{y)G[z)Kl\x-y,z)Df{x-y,z)F{af,Pf,p^) 

J q 

+ g^N, I (z{y)Z{z)K^^{x-y,z)D^^{x,y,z)F{a\\Pl^-p^)+^kf^+eKf^^Z{x) 


-1 


with the new kernel given by 


(28) 


(29) 


k^z’^{x\y,z) = 


_ {x^jy + z) + x'^ (9y^ - Ayz + 9z^) -b x {-9y^ + y'^z + yz'^ - 9z^) - {y - zf {y^ F z^)) 3 g 3 g_ _2 


242,2^2^2 


Z{z){x — z) — 2x{y — 5z) + {y — z)'^) 
24:xy'^ z"^ Z {x) 

Z{y){x - y) -b 2x{py - z) F (y - zf) 
2Axy‘^z‘^Z{x) 


Z{x)F{al\pl^-,p^) 

Z{x)F{al^pl>^-p^). 




(30) 


The renormalization constants have all been dropped in 
favor of the renormalization group improvement terms. 
Note that the renormalization group improvement terms 
for two terms also include an additional factor Z{x). The 
values for the exponents are given in Tab. |llj 

A particular problem of the gluon propagator DSE is the 
appearance of spurious divergences. There are several 
ways to get rid of them. We adopt the one from Ref. [35] 
here. There it was shown that their origin is purely per¬ 
turbative. Consequently we can calculate them analyt¬ 
ically and use these expressions to subtract them. The 
generic structure of the subtraction term is Csuh/p^, with 


Csub = A^CD&^”'”"' V 


n=0 


(ln( AVA|cD))~ ^^" 

n\(—"f + n) 


(31) 


where A is the UV cutoff and Aqcd is defined as the po¬ 
sition of the one-loop Landau pole, viz., Aq^^^ = s 
w is given by 11 Q!(s)/12/7r, where a(s) is the running 

coupling at a perturbative scale s. The coefficient b is de¬ 
termined from the high momentum behavior in the gluon 



9^-17 

3(5-3) 

af 

25 

5-3 


4(35-2) 

3(5-3) 

af 

2(35+13) 

3(5-3) 


2(5+3) 
_5-3 




TABLE II. The exponents a for the RG improvement terms. 
The P’s are all 0 because the exponents were derived as if the 
ghost dressing fnnction were IR constant. 


propagator DSE: 


b = 


6G(s) 


■(G(s)2+2“fz(s)"^“ 

2of ^(^)2-,2/3?^ - 3 CG(s)“3^Z(s)^3* 


)^ 


(32) 


Again the part proportional to does not contribute. 
To deal with the logarithmic divergences a MOM scheme 
is used where the dressings are fixed at the value of the 
Landau gauge result at a perturbative scale. We tested 
explicitly that the combination of these two procedures 
leads to cutoff independent results by varying A^ by a 
factor of 10. 
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Other ingredients for the numerical calculation are func¬ 
tions for extrapolating the dressings beyond the regime 
where they are calculated. For the UV the perturbative 
expressions 


= G{s) 



s 


(33) 




Z{s) 



A2 

A2 

^'^QCD 



(34) 


are used. The anomalous dimensions can be found in 
Tab.|Tj At low momenta we employ a simple extrapolation 
of the form ap^ for the gluon propagator. For the ghost 
propagator we use a constant for ^ = 0 and a/{b + Inp^) 
for ^ > 0. 
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